DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


Shallow  Water  Propagation 

William  L.  Siegmann 
Rensselaer  Polytechnic  Institute 
110  Eighth  Street 

Jonsson-Rowland  Science  Center  1C08 
Troy,  New  York  12180-3590 

phone:  (518)  276-6905  fax:  (518)  276-2825  email:  siegmw@rpi.edu 

Kara  G.  McMahon 
Rensselaer  doctoral  student 

Award  Numbers:  N000140410016 
N000140910638  (Ocean  Acoustics  Graduate  Traineeship) 
http://www.math.rpi.edu/www/ocean_acoustics 


LONG-TERM  GOALS 

Develop  methods  for  deterministic  and  stochastic  acoustic  calculations  in  complex  shallow  water 

environments,  specify  their  capabilities  and  accuracy,  and  apply  them  to  explain  experimental  data  and 

understand  physical  mechanisms  of  propagation. 

OBJECTIVES 

(A)  Treat  propagation  from  narrowband  and  broadband  sources  over  elastic  and  poro-elastic  sediments, 
and  incorporate  realistic  bathymetric,  topographic,  and  geoacoustic  variations. 

(B)  Quantify  acoustic  interactions  with  physical  features  in  the  ocean  volume  and  with  geoacoustic 
features  of  the  ocean  sediment,  and  analyze  and  interpret  experimental  data. 

APPROACH 

(A)  Develop  efficient  and  accurate  parabolic  equation  (PE)  techniques  for  propagation  through 
heterogeneous  sediments.  Treat  range  dependence  and  sediment  layering  by  single  scattering  and 
energy  conservation  methods.  Benchmark  results  using  data  and  special  high-accuracy  solutions. 

(B)  Construct  representations  for  ocean  environmental  and  geoacoustic  variability  using  data  and 
parametric  models.  Determine  acoustic  fields  with  PE,  normal  mode,  and  other  approximation 
methods.  Use  experimental  data  and  computational  results  to  assess  propagation  mechanisms. 

•  Principal  collaborators  are:  Rensselaer  PhD  students,  Dr.  Michael  Collins  (NRL),  Drs.  James 
Lynch,  Timothy  Duda,  and  Ying-Tsong  Lin  (WHOI),  Dr.  Allan  Pierce  (BU,  retired),  and  recent 
Rensselaer  PhD  graduates.  Dr.  William  Carey  (deceased)  was  an  active  collaborator. 
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WORK  COMPLETED 


(A)  Propagation  model  development 

(1 )  Accurate  calculations  range-dependent  sediments  and  applications 

•  High  fidelity  data  from  propagation  over  an  elastic  slab  with  variable  bottom  slopes  in  a  NRL 
tank  experiment  validates  accuracy  [1]  of  a  new  PE  method  designed  for  problems  with  range- 
dependent  bathymetry,  variable  thickness  sediment  layers,  and  topographic  variations  for 
applications  to  beach,  island,  and  coastal  problems. 

•  Additional  verification  of  the  method  is  obtained  by  comparing  results  [2]  between  accuracy 
benchmarks  and  computations  for  environments  with  large  sound  speeds  changes  and  with 
waves  on  range-dependent  elastic  interfaces,  and  guidelines  for  choosing  computational 
parameters  are  provided. 

•  An  application  of  the  new  PE  method  is  made  to  variable  ice  cover  in  polar  regions  [3],  where 
particularly  strong  effects  on  transmission  loss  occur  when  the  surface  ice  terminates  and/or 
regenerates  along  the  propagation  direction. 

•  Another  application  is  to  shallow  regions  with  partially  consolidated  sediment  layers  that  are 
thin  and  have  low  shear  speeds  [4],  which  can  treated  by  modeling  the  layers  as  transitional 
interfaces  and  enforcing  suitable  conditions. 

(2)  New  calculations  for  elastic  and  poro-elastic  sediments 

•  Range-dependent  transversely  isotropic  elastic  sediments,  which  are  a  feature  of  many  coastal 
regions,  are  treated  by  a  PE  new  formulation  [5]  that  enables  evaluation  of  the  relative 
significance  of  anisotropic  effects  on  propagation. 

•  An  initial  propagation  model  that  can  treat  weak  range  dependence  in  transversely  isotropic 
poro-elastic  sediments  demonstrates  that  PE  methods  are  feasible  for  these  environments  [6], 
and  also  indicates  the  influence  of  anisotropy. 

•  Propagation  variables  and  computational  techniques,  developed  for  elastic  sediments  and 
designed  for  range-dependent  layered  environments,  are  generalized  to  multi-layered  poro- 
elastic  sediments  [7],  for  which  computational  results  agree  with  benchmarks  for  range- 
independent  cases. 

(B)  Propagation  mechanism  assessment 

(1)  Nonlinear  internal  wave  effects 

•  If  an  acoustic  mode  propagates  adiabatically  across  a  nonlinear  internal  wave  (NIW)  at  small 
incident  angles  with  the  wave  front,  the  interaction  may  produce  horizontal  Lloyd  mirror 
interference  patterns  [8],  as  predicted  by  other  researchers  and  observed  recently,  and  the 
patterns  have  especially  interesting  features  when  the  NIW  front  has  curvature. 

•  Parameters  that  characterize  NIWs  are  estimated  from  satellite  SAR  images  using  edge 
identification  and  related  techniques  [9]  in  order  to  develop  improved  predictions  of  their 
acoustic  effects,  and  parameter  estimates  are  validated  by  comparisons  with  mooring  and  other 
available  data. 

•  An  adiabatic-mode  transport  theory  is  used  to  develop  a  scattering  model  for  acoustic  energy 
inside  a  NIW  duct  [10]  in  which  wave  front  segments  are  treated  as  scattering  elements,  and  a 
modified  diffusion  equation  describes  how  averaged  intensity  evolves  in  the  duct  and  depends 
on  environmental  and  acoustic  parameters. 
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•  Calculations  from  the  modal  transport  theory  are  used  to  detennine  how  intersecting  NIW 
fronts  influence  average  intensity  [11],  depending  on  the  orientations  of  the  fronts  and  on  the 
asymmetry  of  individual  scattering  elements. 

(2)  Propagation  dependence  on  intrinsic  sediment  attenuation  for  sandy-silt  sediments 

•  A  parametric  description  of  simple  shallow-ocean  waveguides  pennits  development  of  new 
approximations  for  modal  attenuation  coefficients  [12],  and  their  frequency  behavior 
corresponds  to  that  calculated  for  more  complex  waveguides. 

•  Comparisons  between  modal  attenuation  coefficients,  obtained  from  Gulf  of  Mexico  data  and 
from  new  approximations  [13],  were  performed  after  identifying  a  relevant  subset  of  measured 
sound  speed  profiles,  and  the  physical  mechanism  for  the  coefficient  variations  is  explained. 

•  Modal  attenuation  coefficients  and  the  overall  attenuation  of  reduced  transmission  loss  are 
examined  [14]  to  reveal  the  influence  of  sound  speed  gradients  at  different  locations  in  the 
water  column  and  of  sound  speed  profiles  in  the  sediment. 

•  Relatively  compact  formulas  for  averaged  transmission  loss  in  range-independent  waveguides 
are  derived  from  mode  theory  [15],  and  they  reduce  to  well-known  results  of  Rogers  and  others 
for  high  frequencies  and  either  isospeed  or  constant-gradient  water  sound  speed  profiles. 

(3)  Card-house  structure  theory  for  muddy  sediments 

•  One  feature  of  this  model  is  the  presence  of  electric  charges  on  bubbles,  and  under  unifonn 
conditions  it  is  found  that  charged  bubbles  are  required  to  be  non-spherical  [16],  which  has 
been  observed  in  laboratory  experiments  and  is  anticipated  for  mud  in  situ. 

•  Another  key  model  component  is  that  based  on  their  chemical,  electrical,  and  material 
structure,  mud  platelets  are  hypothesized  to  behave  like  electric  quadrupoles,  so  that  estimates 
for  shear  wave  speed  are  obtained  [17]  by  computing  the  oscillation  frequency  of  a  hinged  joint 
fonned  by  two  platelets  which  interact  end  to  side. 

•  The  shear  speed  calculation  is  extended  to  account  for  platelets  for  which  the  interaction  is 
modeled  as  a  cantilevered  beam  rather  than  a  hinged  joint  [18],  and  estimates  are  consistent 
with  low  shear  speed  values  found  in  high-porosity  mud  sediments. 

•  Estimates  of  another  important  mud  parameter,  porosity,  are  found  by  hypothesizing  that  mud 
platelets  aggregate  by  an  idealized  process  of  fractal  type  [19],  which  is  validated  by  finding  a 
fractal  dimension  for  the  process  and  comparing  with  a  result  detennined  by  an  independent 
computation. 

RESULTS  (from  two  selected  investigations) 

(A)  A  critical  capability  for  ocean  acoustic  data  analysis  and  other  applications  is  efficient  and  accurate 
propagation  calculations  in  shallow  water  waveguides  with  range-dependent  poro-elastic  sediment 
layers.  Porosity  and  elasticity  are  important  physical  mechanisms,  because  energy  transfers 
between  compressional  and  shear  modes  produce  significant  acoustic  intensity  and  phase  changes, 
and  because  attenuation  may  also  increase  substantially.  A  major  computational  challenge  arises 
because  wave  number-energy  spectra  are  much  broader  for  poro-elastic  sediments  than  for  fluids. 
Our  approach  is  based  on  the  Biot  poro-elastic  theory,  which  is  the  most  common  model  for  ocean 
sediments.  The  only  published  poro-elastic  PE  appeared  over  a  decade  and  a  half  ago,  and  even 
with  an  extension  for  sediment  anisotropy  [5],  its  computational  capabilities  are  too  limited  for 
applications.  It  is  essential  that  recent  progress  for  elastic  sediments  be  generalized  to  layered 
poro-elastic  media.  For  example,  one  of  the  new  methods  has  been  validated  with  high  quality  data 
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from  an  NRL  experimental  series  using  elastic  slabs  [1]  and  with  accurate  computations  for 
environments  with  range  dependent  interface  waves  [2].  These  methods  were  developed  from 
advances  including  formulations  with  different  dependent  variables,  rotations  of  coordinates  at 
range  locations  of  bathymetry  slope  changes,  and  improved  range-dependent  corrections  at  stair¬ 
step  approximations  for  ocean  and  sediment  parameters.  Last  year  we  reported  progress  for  the 
poro-elastic  PE  by  generalizing  its  original  variables  to  another  set  (denoted  by  q  ),  two  of  which 
appear  in  new  elastic  results.  The  top  left  panel  in  Figure  1  illustrates  that  loss  calculations  for 
poro-elastic  sediments  with  one  or  two  layers  can  be  drastically  different.  We  found  that  the 
method  is  accurate  for  one-layer  sediments;  results  are  not  displayed  here.  However,  it  is  possible 
for  accuracy  to  deteriorate  unacceptably  for  two-layer  sediments,  as  shown  in  the  top  right  panel. 
This  dilemma  prompted  a  full  reexamination  of  formulations  for  poro-elastic  media,  with  the  result 
that  a  new  set  of  variables  was  devised  for  layered  and  range-dependent  environments  [7].  The 
middle  panel  shows  calculations  for  a  range-independent  environment  with  the  revised  formulation 
(denoted  by  q  ),  along  with  a  benchmark  solution  from  the  wavenumber  integration  code  OASIS. 
The  two  curves  are  in  excellent  agreement.  This  example  and  others  demonstrate  that  the  q 
formulation  produces  accurate  TL  results  for  multi-layered  poro-elastic  media.  Moreover,  the 
formulation  has  capabilities  which  are  needed  for  applications.  The  bottom  left  panel  shows  loss 
color  contours  for  a  source  close  to  an  interface  between  water  and  a  poro-elastic  sediment.  A 
distinctive  interference  pattern  produced  by  an  interface  wave  and  a  water  mode  is  apparent.  The 
Scholte-like  interface  wave  is  directly  analogous  to  the  important  Scholte  wave  at  a  water-elastic 
interface.  This  is  shown  by  the  bottom  right  panel,  where  a  horizontal  wave  number  spectrum 
indicates  that  the  spectral  position  of  the  interface  wave  is  near  that  of  a  true  Scholte  wave  when 
the  poro-elastic  sediment  becomes  elastic.  From  these  and  other  calculations,  we  conclude  that  the 
q  PE  solution  has  new  capabilities  for  accurate  propagation  calculations  in  multi-layered  poro- 
elastic  environments.  Future  work  includes  improving  capabilities,  accuracy  and  efficiency  for 
range-dependent  problems,  and  additional  benchmarking. 


4 


Figure  1  (next  page).  Accurate  calculation  of  propagation  in  layered poro-elastic  sediments  is 
obtained  from  an  appropriate  formulation  of physical  variables  in  the  PE  method.  Two  range- 
independent  environments  are  described  here,  both  having  a  100  m  isospeed  ocean  with  cw  =  1500 
m/s  overlying  a  poro-elastic  sediment.  Environment  A  has  a  single  deep  layer  with  geoacoustic 
parameters:  density  solid  ps  =  1.3  gm/cni  ' ,  porosity  a  =  0.2,  sound  speeds  cpi  =  1 700  m/s,  cP2  =  1250 
m/s,  and  cs  =  800  m/s,  and  attenuations  ppl  =  0.1  dB/) ,  pp2  =  1  dB/) ,  and /f  =  0.2  dB/) .  Environment 
B  has  a  20  m  thick  upper  sediment  layer  that  is  the  same  as  in  Environment  A  and  overlies  a  deep 
layer  with  geoacoustic paramters:  ps  =  1.7 gm/cni ' ,  a  =  0.2,  cpl  =  2040  m/s,  cp2  =  1250  m/s,  cs  =  950 
m/s,  ppi  =  0.3  dB//. ,  pp2  =  5  dB/). ,  and  ps  =0.5  dB /) .  Higher  values  of  pp2  reflect  higher  attenuations 
for  the  Biot  slow  compressional  wave.  In  Figs.  l(a)-(c),  transmission  loss  (TL)  curves  (dB  re:  1  m) 
are  shown  between  30  and  90  dB  along  a  2  km  track  which  is  75  m  below  the  upper  surface,  for  a  50 

Hz  source  at  75  m  depth. 

Fig.  1(a)  shows  curves  calculated  by  an  earlier  PE  formulation  with  variables  q  =  (ux,  w,  C ),  which 
generalize  a  successful  elastic-media  method.  The  dashed  TL  curve  for  Environment  A  shows  a 
two-mode  pattern,  while  the  solid  curve  for  Environment  B  shows  one  or  more  additional  modes  that 
produce  a  shorter  pattern  wavelength  of  about  300  m.  The  dissimilar  TL  curves  mean  that  a  single 
poro-elastic  layer  does  not  accurately  represent  an  environment  with  two  or  more  layers,  when  the 
upper  layer  is  one  wavelength  thick.  In  Fig.  1(b)  the  solid  TL  curve  is  the  same  as  in  Fig.  1(a).  The 

dashed  curve  also  uses  Environment  B  and  is  calculated  from  the  wave  number  integration 
implementation  OASIS.  The  curves  have  similar  shapes  but  differ  in  amplitude  and  pattern  phase. 

Fig.  1(c)  curves  use  the  two-layer  Environment  B.  They  are  plotted  over  a  5  km  track  and  evolve  to  a 
three-mode  pattern.  The  solid  curve  is  calculated  by  a  formulation  with  a  new  set  of  dependent 
variables  q  =  (T ,  w,  A  ),  and  the  dashed  curve  is  the  solution  from  OASIS.  The  two  curves  are  in 
excellent  agreement  and  demonstrate  that  the  new  formulation  produces  accurate  TL  results  for 

multi-layered  poro-elastic  media. 

Fig.  1(d)  illustrates  that  the  new  formulation  has  useful  capabilities,  such  as  calculating  the 
influence  of  interface  waves  on  TL.  The  same  two-layer  sediment  Environment  B  is  used  with  a  25 
Hz  source  at  98  m,just  above  the  interface.  Color  TL  contours  over  2  km  in  range  and  between  70 
and  150  m  depth  are  shown  with  a  dynamic  range  of  40  dB.  In  addition  to  the  mainly  three-mode 
pattern,  the  contours  show  a  rapid  oscillation  above  and  below  100  m  from  an  interface  wave. 

Fig.  1(e)  provides  the  horizontal  wave  number  spectrum  between  0.05  and  0.25  /m,  normalized  with 
maximum-peak  amplitude  one,  for  the  calculation  in  Fig.  1(d)  at  a  receiver  depth  98  m.  One  large 
and  two  smaller  peaks  between  0.07  and  0.11  /m  correspond  to  the  three  water  modes,  and  another 
peak  (with  half  the  maximum  amplitude)  near  0.23  /m  corresponds  to  the  poro-elastic  analog  of  a 
Scholte  interface  wave.  This  Scholte-like  wave  has  a  speed  near  695  m/s,  which  is  about  5%  larger 
than  the  Scholte  wave  speed  if  the  sediment  is  treated  as  elastic.  From  these  and  other  calculations, 
we  conclude  that  the  q  solution  has  the  capability  for  accurate  propagation  calculations  in  multi¬ 
layered  poro-elastic  environments. 
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(B)  Both  scientific  understanding  and  specific  knowledge  of  acoustic  properties  of  shallow  water 
bottoms  are  fundamental  for  many  applications  of  interest.  The  widespread  sandy-silt  bottoms 
have  received  most  attention  and  are  handled  effectively  by  elastic  or  Biot  poro-elastic  models. 
However,  high-porosity  seabed  mud  is  different,  because  neither  of  these  models  correctly  predicts 
or  explains  key  geoacoustic  properties  of  porosity,  shear  and  compressional  sound  speeds,  and 
attenuations.  A  card-house  structure,  recently  proposed  by  Allan  Pierce  and  William  Carey,  can 
potentially  revolutionize  approaches  for  treating  mud  bottoms.  The  structure  relies  on  the  fact  that 
mineral  components,  which  are  in  the  form  of  thin  platelets,  have  a  particular  distribution  of 
electric  charges  that  causes  them  to  resemble  small  quadrupole  sheets.  This  leads  to  repulsion  for 
end-to-end  or  face-to-face  contact  and  attraction  between  ends  and  faces,  which  in  turn  causes 
platelets  to  aggregate  into  structures  of  card-house  fonn.  The  aggregates  may  have  pores  where 
bubbles  of  non-spherical  shape  can  occur,  and  an  initial  explanation  for  such  bubble  shapes  has 
been  proposed  [16].  In  addition,  different  models  for  the  electrical  and  mechanical  interactions  of 
platelets  provide  mechanisms  for  estimating  shear  speeds,  depending  on  whether  the  platelets  are 
predominantly  rigid  [17]  or  elastic  [18].  Developments  along  these  lines  are  critical  to  demonstrate 
the  effectiveness  of  the  card-house  model  for  applications.  In  addition,  the  evolution  and 
geometric  properties  of  structures,  and  estimation  of  volume  fractions  of  platelets  and  pores 
(porosity  [1),  are  themselves  important  to  resolve.  We  considered  known  aggregation  processes  and 
concluded  that  for  mud,  the  most  likely  is  cluster-cluster  aggregation  constrained  mainly  by  platelet 
diffusion.  A  collection  of  physically-based  rules  are  hypothesized  for  such  a  process,  along  with 
several  computational  simplifications  such  as  two-dimensional  platelet  interactions  and  horizontal 
or  vertical  platelet  orientations.  The  top  left  panel  of  Figure  2  illustrates  a  random  initial 
distribution  of  platelets  on  a  unifonn  rectangular  grid  with  point  spacing  so  that  a  platelet  fits  on 
three  points.  Computational  iterations  are  performed  that  allow  randomly  selected  platelets  to 
translate  or  rotate.  The  top  right  panel,  after  1000  iterations,  shows  clusters  that  have  formed  and 
suggests  their  increasing  complexity  as  they  interact.  The  lower  left  panel  validates  that  the 
aggregation  process  has  a  fractal  nature.  The  fractal  dimension  D  should  be  the  slope  in  a  log-log 
plot  of  the  number  of  platelets  in  a  cluster  versus  an  effective  (dimensionless)  cluster  size, 
calculated  from  a  radius  of  gyration.  The  slope  D  =  1 .60  is  significantly  less  than  two  which  is  the 
value  for  non-fractal  behavior.  The  lower  right  panel  involves  the  fractional  solid  content  1  -  P  of 
a  cluster,  calculated  from  a  ratio  of  estimated  solid  and  total  areas.  This  log-log  plot  uses  a 
specific  combination  as  ordinate  and  the  radius  as  abscissa  in  order  that  its  slope  is  another  estimate 
for  D  =  1.56.  The  closeness  of  the  two  D  values  is  additional  validation  that  the  aggregation 
process  is  fractal.  Moreover,  if  a  selection  is  made  for  platelet  aspect  ratio,  such  as  0.05  for  the 
mineral  kaolinite,  an  estimate  for  porosity  P  is  in  the  range  0.89  -  0.94.  This  is  consistent  with 
accepted  values  near  0.90  for  mud  for  which  the  card -house  theory  should  apply.  We  conclude 
from  these  and  other  computations  with  the  preliminary  model  that  a  fractal  hypothesis  for  card- 
house  aggregation  is  reasonable,  and  that  estimates  for  fractal  dimension  are  consistent  with 
experimental  observations.  More  importantly  for  geoacoustics,  porosity  estimates  obtained  from 
the  computations  are  consistent  with  experimental  observations.  Future  work  includes  relaxing 
constraints  in  the  computational  model,  continuing  the  determination  of  porosity  and  sound  speed, 
and  specifying  sensitivities  of  parameter  estimates. 
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Figure  2  (next  vase).  The  card-house  theory  for  mud  structure  should  permit  estimates  for  key 
geoacoustic  parameters  such  as  porosity  and  shear  sound  speed.  A  computational  approach  based 
on  modeling  aggregation  of  platelets  provides  preliminary  estimates  of  porosity.  The  mineral 
constituents  of  mud  are  modeled  as  thin  platelets  with  length  L,  thickness  h,  and  small  aspect  ratio  S 
=  h/L.  Electrically  they  resemble  quadrupoles  that  repel  face-to-face  and  end-to-end  but  attract 
face-to-end.  Aggregation  is  treated  as  a  diffusion-limited  cluster-cluster  process  that  follows 
certain  modeling  rules,  principally  (1)  interactions  are  two-dimensional,  (2)  platelets  always  lie 
horizontally  or  vertically  on  a  uniform  grid  of  points  with  spacing  equal  to  L/2,  (3)  platelets  move 
horizontally,  vertically,  or  by  rotating  about  their  centers,  (4)  end-to-end  and  face-to-face  platelet 
contact  and  overlapping  platelets  are  forbidden,  and  (5)  platelets  undergoing  center-to-end  contact 

remain  attached  permanently. 

Fig.  2(a)  is  an  example  of  a  random  initial  placement  of  nearly  1500  platelets,  shown  on  a  40  by  40 
grid  with  spacing  0.5  (all  lengths  in  units  of  L).  At  this  stage  very  few  platelets  have  clustered.  A 
series  of  iterations  are  performed,  allowing  randomly  selected  platelets  to  move  up  or  down,  or  left 
or  right,  one  grid  spacing  if feasible.  Other  randomly  selected  platelets  perform  a  rotation  of  90 
degrees  if  feasible.  Fig.  2(b)  shows  the  result  after  1000  iterations.  A  total  of  17  clusters  formed 
with  a  variety  of  sizes,  several  of  which  escaped  the  area  shown.  In  the  current  version  of  the 
process,  the  clusters  typically  have  pores  that  are  accessible  to  the  grid. 

Fig.  2(c)  is  a  demonstration  of  the  fractal  nature  of  the  hypothesized  aggregation  process.  A 
number  of  experimental  investigations  concluded  that  mud  aggregates  possess  fractal  geometry, 
with  Minkowski-Bouligand  (box  counting)  dimension  of  1.5  or  somewhat  higher.  A  reasonable 
expectation  is  that  a  two-dimensional  non-fractal  cluster  has  dimensionless  mass  M  satisfying 
M  oc  Rg2  where  Rg  is  a  calculated  radius  of  gyration  divided  by  L.  Because  the  number  N  of 
platelets  in  a  cluster  is  proportional  to  the  dimensionless  mass,  N  oc  R2.  These  relations  generalize 
to  M  oc  Rg°  and  N  oc  Rg°  for  a  fractal  cluster,  where  D  is  the  fractal  dimension.  The  figure  shows  a 
plot  of  log  N  (from  about  2.19  to  3. 05)  vs.  log  Rg  (from  about  0. 69  to  1.20)  for  approximately  three 
dozen  aggregates,  all  of  which  were  obtained  from  model  computations  and  contain  at  least  150 
platelets.  The  best  fit  line  (solid  red)  has  slope  1.60,  which  is  an  estimate  for  D  that  is  consistent 
with  previous  work.  A  comparison  line  for  non-fractal  behavior  (dashed  blue)  has  slope  2. 

Fig.  2(d)  provides  results  for  aggregate  porosity  (3.  The  fractional  solid  content  of  a  cluster,  taken 

n _ o 

as  proportional  to  its  (dimensionless)  density,  is  expected  to  satisfy  the  relation  1  —  (3  oc  Rg 
The  left  side  can  be  calculated  by  estimating  solid  and  total  areas  for  a  cluster  using  a  box  counting 
procedure.  It  is  worth  noting  from  this  relation  that  S  appears  only  in  the  proportionality  factor, 
and  consequently  the  exponent  D  is  independent  of  S.  The  figure  shows  a  plot  of  log  ((1  -  (3)Rg  ) 
(from  about  0.37  to  1.22)  vs.  log  Rg  (from  about  0. 69  to  1.20),  which  has  a  slope  D  =  1. 56.  The 
closeness  of  the  two  estimates  for  fractal  dimension  is  important  for  the  validity  of  the  original 
hypothesis.  With  D  given,  porosity  may  be  obtained  by  selecting  a  value  for  S.  For  kaolinite 
platelets  a  typical  value  is  S  =  0.05,  producing  porosity  estimates  in  the  range  (3  =  0.89-  0.94.  These 
are  consistent  with  accepted  values  near  0.90  in  mud  for  which  the  card-house  theory  should  apply. 
We  conclude  from  computations  with  this  initial  model  version  that  a  fractal  hypothesis  for  card- 
house  aggregation  is  reasonable  and  that  estimates  for  fractal  dimension  are  consistent  with 
experimental  observations.  More  importantly  for  geoacoustics,  porosity  estimates  obtained  from  the 

model  are  consistent  with  experimental  observations. 
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IMPACT/APPLICATIONS 


New  or  enhanced  capabilities  are  provided  for  propagation  predictions  that  depend  on  physical 
properties  of  shallow  water  sediments,  including  layering,  elasticity,  porosity,  and  anisotropy.  Range- 
dependent  variability  from  bathymetry,  topography,  and  sediment  layer  interfaces  in  propagation 
calculations  can  be  treated.  Intensity  attenuation  and  coherence  statistics  that  result  from 
environmental  fluctuations  and  other  experimental  variability  can  be  found  more  efficiently.  Data 
analyses  and  model  comparisons  allow  specification  for  application  purposes  of  the  relative 
significance  of  key  physical  mechanisms:  linear  versus  nonlinear  frequency  dependence  of  sediment 
attenuation,  sediment  heterogeneity  versus  homogeneity,  water  column  versus  bathymetric  variability, 
water  column  scattering  versus  refraction,  and  vertical  versus  horizontal  mode  coupling  due  to 
nonlinear  internal  waves  and  bathymetry.  Results  from  modeling  and  data  analyses  of  experiments, 
particularly  experiments  off  the  New  Jersey  Shelf,  are  partly  aimed  toward  improving  shallow  water 
sonar  systems  and  predictions.  Propagation  model  implementations,  analysis  tools,  and  data 
representation  techniques  are  distributed  to  university,  laboratory,  and  research/development  groups. 

RELATED  PROJECTS 

•  Continuing  projects  with  Dr.  Michael  Collins  include  [5]  and  a  monograph  on  new  parabolic  wave 
equation  models  and  applications  [21],  for  which  the  principal  research  results  are  nearly  complete 
and  chapter  drafts  are  prepared. 

•  In  addition  to  investigations  with  Drs.  James  Lynch,  Timothy  Duda,  and  Y.-T.  Lin  and  their 
colleagues  [9]-[  1 1]  of  propagation  effects  from  waveguides  of  variable  structure  generated  by 
nonlinear  internal  waves,  a  related  project  characterizes  whispering  gallery  and  other  modes  [20]. 

•  Research  with  Dr.  Allan  Pierce  [15]-[19]  focuses  on  propagation  variability  from  sediment 
geoacoustic  structure  and  attenuation  and  on  quantifying  predictions  from  his  card-house  theory  of 
mud  structure  and  acoustics. 
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